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Motivated by the experimentally observed shear-induced destabilization and reorientation of smec- 
tic A like systems, we consider an extended formulation of smectic A hydrodynamics. We include 
both, the smectic layering (via the layer displacement it and the layer normal p) and the director 
n of the underlying nematic order in our macroscopic hydrodynamic description and allow both 
directions to differ in non equilibrium situations. In an homeotropically aligned sample the nematic 
director does couple to an applied simple shear, whereas the smectic layering stays unchanged. This 
difference leads to a finite (but usually small) angle between h and p, which we find to be equivalent 
to an effective dilatation of the layers. This effective dilatation leads, above a certain threshold, to 
an undulation instability of the layers. We generalize our earlier approach [Rheol. Acta 39 (3), 15] 
and include the cross couplings with the velocity field and the order parameters for orientational 
and positional order and show how the order parameters interact with the undulation instability. 
We explore the influence of various material parameters on the instability. Comparing our results 
to recent experiments and molecular dynamic simulations, we find a good qualitative agreement. 

PACS numbers: 61.30-v, 47.20Ft, 83.50Ax, 05.70Ln 
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I. INTRODUCTION 



Submitted to an applied shear flow, many complex liq- 
uids show an interesting coupling between their internal 
structure and the flow field. For smectic A like systems 
(including block copolymers, lyotropic systems and side- 
chain liquid crystalline polymers) this coupling may in- 
duce reorientation of the layers. Experiments on a variety 
of systems which differ significantly in their microscopic 
details show nevertheless striking similarities in their 
macroscopic behavior under shear. The systems under 
investigation include block copolymers @- Jli Hj^PIj i|i 
low molecular weight (LMW) liquid crystal s |H H ff, ly- 
otropic lamellar phases (both LMW [0 |ll|, [l^and poly- 
meri c [|l3| ) , and liquid crystalline side-chain polymers 
p^ , |l5[. These experiments use either a steady shear 
(typically for the low viscosity systems e.g. in a Couette 
cell) or large amplitude oscillatory shear (often in the 
highly viscous polymeric systems, e.g. in a cone-plate or 
plate-plate geometry). Due to these experimental differ- 
ences a direct comparison between the different systems 
is not always straightforward. The common features of 
all these experiments can be described as follows. Start- 
ing with an aligned sample where the layers are parallel 
to the planes of constant velocity ( "parallel" orientation) , 
the layering is stable up to a certain critical shear rate 
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At higher shear rates two different 
situations are observed. Depending on the system either 
multi-lamellar vesicles ||l^, |l^, ^ ("onions", typically 
in lyotropic systems) or layers perpendicular to the vor- 
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ticity direction § |, |, |, |, |, |, Q ("perpendicular" 
orientation, typically in solvent free systems) form. In 
some of the systems a third regime is observed at even 
higher shear rates with a parallel orientation |lOj. If 
the starting point is rather a randomly distributed lamel- 
lar phase, the first regime is not observed ||l|, |[ |l^, p^ . 
This last point illustrates that experiments on layered liq- 
uids depend on the history of the sample. In our further 
discussion we will restrict ourselves to systems showing a 
well aligned parallel orientation before shear is applied. 

The experimental similarities between different sys- 
tems indicate, that the theoretical description of these 
reorientations can be constructed, at least to some ex- 
tent, from a common basis independent of the actual 
system (on the other hand, a description including the 
differences between the systems under investigation must 
refer closer to their microscopic details). When looking 
for a macroscopic description, the well established stan- 
dard smectic A hydrodynamics ^ ^ ^ is a good 
starting point for such a theoretical approach. 

As first shown by Delaye et al. and Clark and 

Meyer thermotropic smectic A liquid crystals are 
very sensitive against dilatations of the layers. Above a 
very small, but finite, critical dilatation the liquid crystal 
develops undulations of the layers to reduce the strain 
locally. Oswald and Ben-Abraham considered dilated 
smectic A under shear . When a shear flow is applied 
(with a parallel orientation of the layers), the onset for 
undulations is unchanged only if the wave vector of the 
undulations points in the vorticity direction (a similar sit- 
uation was later considered in |Q ) . Whenever this wave 
vector has a component in the flow direction, the onset 
of the undulation instability is augmented by a portion 
proportional to the applied shear rate. No destabilizing 
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mechanism for well aligned parallel layers is present in 
the standard smectic A hydrodynamics. 

Recently we proposed an extended hydrodynamic de- 
scription of smectic A liquid crystals. Using both, 
the director of the underlying nematic order and the layer 
normal of the smectic layers, we showed the possibility 
of a shear-induced undulation instability in well aligned 
parallel layers. Within the framework of irreversible ther- 
modynamics (which allows the inclusion of dissipative as 
well as reversible effects) we derived macroscopic hydro- 
dynamic equations for the system and performed a lin- 
ear stability analysis of these equations (using a number 
of approximations). As always, a linear stability analy- 
sis is limited to the onset of the first instability. Other 
theoretical approaches to these reorientation phenomena 
have been undertaken by Bruinsma and Rabin p7| , Zil- 
man and Granek p^ ] (both papers are considering the 
influence of the shear on layer fluctuations) and Williams 
and MacKintosh [ p9[ (minimizing a free energy density 
including a coupling to the applied shear stress). To our 
knowledge, no macroscopic hydrodynamic approach be- 
sides |6| has been published up to now. 

The present paper is structur ed a s follows: After a 
brief review of the model in Sec. II A and its implemen- 
tation in Sec. II B we extend the basic model of ES, Efll 



in the following sections. Especially we include the cross 
coupling to the velocity field and the moduli of the ne- 
matic and smectic order parameters. It turns out that 
the coupling terms to the velocity are important since 
they can change the critical parameters significantly. We 
find that the moduli of the order parameters also show 
undulations and, thus, regions with a reduced order para- 
meter can be identified. The comparison of the different 
levels of approximations shows that the basic model is 
contained in this more general analysis as a special case. 
We also compare our results to experiments and molec- 
ular dynamic simulations and show that an oscillatory 
instability is extremely unlikely to occur. 



II. MODEL AND TECHNIQUE 
A. Physical idea of the model 

In a smectic A liquid crystal one can easily define two 
directions: the normal to the layers p and an average 
over the molecular axes, the director, h. In the stan- 
dard formulation of smectic A hydrodynamics these two 
directions are parallel by construction. Only in the vicin- 
ity of phase transitions (either the nematic-smcctic A or 
smectic A-smectic C*) it has been shown that director 
fluctuations are of physical interest |30| Nev- 
ertheless n and p differ significantly in their interaction 
with an applied shear flow. 

We consider a situation as show in Fig. |^. A well 
aligned smectic A liquid crystal is placed between two 
parallel and laterally infinite plates. The upper plate 
(located at 2; = d/2) moves with a constant velocity 



FIG. 1: At the level of the approximation we use in this paper, 
all experimental shear geometries are equivalent to a simple 
steady shear. We choose our system of coordinates such that 
the normal to the plates points along the z-axis and the plates 
are located at z = 
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FIG. 2: A finite angle Q between n and p leads to a tendency 
of the layers to reduce their thickness. Supposing a constant 
number of layers in the sample, this tendency is equivalent to 
an effective dilatation of the layers. For small angles between 
h and p the relative effective dilation is given by fP' 12. 




Vu = djix/'i to the right and the lower plate (at —d/2) 
moves with the same velocity in the opposite direction 
{vi — —d'jex/2). Thus the sample is submitted to an 
average shear given by (w„ — vi)/d = 'j. As mentioned 
above, a three dimensional stack of parallel fiuid layers 
cannot couple directly to an applied shear flow. Neither 
does the layer normal: it stays unchanged as long as the 
flow direction lies within the layers. In contrast, it is 
well known from nematodynamics that the director ex- 
periences a torque in a shear flow. This torque leads — 
in the simplest case — to a fiow aligning behavior of the 
director. The key assumption in the model of is that 
this torque is still present in a smectic A liquid crystal 
and acts only on the director n and not on the layer nor- 
mal p. An energetic coupling between h and p ensures 
that both directions are parallel in equilibrium. 

Submitted to a shear flow the layer normal will stay un- 
changed, but the director will tilt in the direction of the 
flow until the torques due to the flow and due to the cou- 
pling to the layer normal balance one another. For any 
given shear rate a flnite, but usually small, angle 6 be- 
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FIG. 3: Above a certain threshold the effective dilatation due 
to the director tilt will lead to buckling of the layers. Note 
the difference in directions: the director is tilted in the ffow 
direction, whereas the wave vector points along the y-axis. 
This configuration cancels the direct coupling between the 
flow and the buckling. 
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tween n and p will result. This finite angle has important 
consequences for the layers: Since the preferred thickness 
of the layers is proportional to the projection of the di- 
rector on the layer normal, a finite angle between those 
two directions is equivalent to an effective dilatation of 
approximately 0^/2 (see Fig. H). If we assume a constant 
total sample thickness and exclude effects of defects, the 
system can accommodate this constraint by layer rota- 
tions. A global rotation of the layers is not possible, but 
they can rotate locally (as in the case of dilated smectic A 
liquid crystals ||2l|, ) . This local rotation of the layers 
leads to undulations as shown in Fig. ^. These undu- 
lation are a compromise between the effective dilatation 
(which is not favorable for the system) and the curvature 
of the layers due to the undulations (which costs energy) . 
In the static case of dilated smectic A liquid crystals no 
direction is preferred, but Oswald and Ben- Abraham 
have shown that this symmetry is broken if an additional 
shear is applied to the system. In this case the standard 
formulation of smectic A hydrodynamics predicts that 
the wave vector of the undulations will point along the 
neutral direction of the shear. In this paper we will as- 
sume that this result of Oswald and Ben- Abraham also 
holds in the case of our extended formulation of smectic 
A hydrodynamics (see Fig. |^). 



Let us briefly review the essential ingredients to this 
procedure (for more details of the method see and for 
our model [p5[). For a given system the hydrodynamic 
variables can be split up in two categories: variables re- 
flecting conserved quantities (e.g. the linear momentum 
density, the mass density etc.) and variables due to spon- 
taneously broken continuous symmetries (e.g. the ne- 
matic director or the layer displacements of the smectic 
layers). Additionally, in some cases non- hydrodynamic 
variables (e.g. the strength of the order parameter [^) 
can show slow dynamics which can be described within 
this framework (see e.g. [|o[ |4|). 

Using all these variables the relations, which form the 
starting point for the further calculations, can be con- 
structed. These relations are: the energy density e, the 
dissipation function i?, the Gibbs-, and the Gibbs-Duhem 
relation. To illustrate the idea of our model we split up e 
and R into several parts according to the different origin 
of the variables. 

- conserved quantities (index cons) 

- symmetry variables (index sym) 

- the modulus of the order parameter (index ord) 

In the spirit of our model two order parameters play a 
role: the nematic tensorial order parameter and the 
smectic A complex order parameter $. For practical rea- 
son we use the director n and the modulus S*'"' in the 
uniaxial nematic case [Qij = ^S'^'^^riinj — \5ij)] and the 
layer displacement u and the modulus S^""^ in the smectic 
A case [(p — S'^^^ eyiY>{iqQ{z — u)}]. Here, as in the rest of 
the paper, we refer to the system of coordinates defined 
in Sect. II A. We note that u is only a good variable if 
we consider small deformations of the layers. For large 
layer deformations the phase Lp — z — u is the appropri- 
ate variable |3^, |3^. In our further discussion, we will 
concentrate on the parts due to symmetry variables and 
the order parameters, while for terms already present in 
the isotropic fluid see e.g. [|o[ 

Let us first consider the energy density. The conven- 
tions of notation introduced by the following equations 
are summarized in Tab. ffl. 
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ord 
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B. Implementation of the model 

To generate the macroscopic hydrodynamic equations 
we follow the procedure given by the framework of irre- 
versible thermodynamics |3^. This method has success- 
fully been applied in many cases to derive the macro- 
scopic hydrodynamic equations of complex fluids (see e.g. 
[^,20, Q). The advantage of this method is its sys- 
tematic way of deducing the governing equations. Once 
the set of variables is given, the macroscopic hydrody- 
namic equations follow by applying basic symmetry ar- 
guments and thermodynamic considerations. 



Econs, which is identical to the isotropic fluid, is discussed 
elsewhere |^ , |3^ . The symmetry part reads: 

1 9 

+ -K3 [fi X (V X h)]- 



'-Bo 
2 ° 



V,M+(l-n,)-i(Viu)' 



^Bi {h X p)^ 



(2) 
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symbol explicit form comment 



K 

Bo 



- (n.a) 



Tin) 



K 

Bo 

r (n,s) 



(n) 



niUj 



bending modulus of layers 

compressibility of layers 

coupling between the director 
and the layer normal 

variations of the order 
parameter (nematic and 
smectic, respectively) 

gradients terms of the order 
parameter (nematic) 

cross-coupling between the 
director and order parameter 
(nematic) 

gradients terms of the order 
parameter (smectic) 



TABLE I: Summary of the notation. In these definitions we 
use the transverse Kronecker symbol S:^ — Sij — niUj. Due 
to the thermodynamic stability of the systems the following 
combinations of constants must be positive: Bo, Bi, K, Lo, 



L^i"'"', and - KL\\ 
equivalence of K and Ki. 



For the last relation we used the 



In Eq. (||) the spirit of the model becomes clear. We 
combine the properties of a nematic Uquid crystal (the 
first two lines) with these of a smectic A (the third and 
fourth line) and couple both parts (the last line) in such 
a way that fi and p are parallel in equilibrium. As al- 



ready discussed earlier 



simplifies considerably 



by dropping higher order terms and assuming a small 
angle between h and p. Splay deformations of the di- 
rector are generally considered as higher order correc- 
tions to dilatations of the smectic layers. Twist defor- 
mations arc forbidden in standard smectic A hydrody- 
namics and must be small as long as the angle between 
n and p is small. Additionally, the difference between 
the splay deformation of the director field Ki/2 (V • n)^ 



and bending of the layers K/2 (Vj 



is negligible. 



Consequently we combine splay and bend in a single 
term with a single elastic constant which we call K': 
Ki/2 {V-nf+K/2 {V\uf kK'/2 (V»^ In the 
following we drop the prime and call the new elastic con- 
stant K. The approximated version of Csym is now given 
by: 



^sym 



-K Nlu) 



'-Bo 
2 ° 



-t- (1 - n,) - i (Viu)' 



{fi X pf 



name 


variable 


conjugate 


mass density 


P 




momentum density 
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_^ 
V 


nematic director 


h 


h 


smectic layer displacement 


u 




variation of the modulus of "j 






the order parameter (either > 







nematic or smectic) 



TABLE II: Variables and their conjugates, i.e. the corre- 
sponding thermodynamic forces 



In our model the moduli of the nematic and smectic order 
parameter play similar roles, so we will deal with both. 
Since we consider a situation beyond the phase transition 
regime, the equilibrium value of the order parameter is 
non-zero {Sq"''^\ for both nematic and smectic) and only 
its variations s^"'*) can enter the energy density (S'^"''*-' = 



S, 



{n,s) 




Jn,s) 



An) 
ord 



(4) 



4:^. = \L,{s^^^)\\L[%{v.s^^^){v,si^^) (5) 



By a similar construction we write down the dissipa- 
tion function as (see Tab. |l| for a list of the thermody- 
namic variables and their conjugates) 



R 



Rs 
R. 



ord 



^cons ~t~ Rsyin Rord (6) 

= ^i^ijki (VjWj) (VjUfe) + -Ro (7) 

= ^M.^/i,+Ap*2 (8) 

= ia(^) S(^)' (9) 



where Rq summarizes further terms due to conservations 
laws, which do not enter our further calculation, and (af- 
ter M) 



Vijkl 



+ 2{vi + V2- 2iy3)ntnjnkni 
+ [vs - V2){njni5ik + UjnkSu 

+ - V2)5ij5ki 

+ - + '^2)iSijnkni + SkiniHj) (10) 



As mentioned in Sec. 11 A we consider a shear induced 



(3) smectic C like situation (but with a small tilt angle, i.e 



a weak biaxiality) . We neglect this weak biaxiality in the 
viscosity tensor and use it in the uniaxial formulation 
given above (with the director n as the preferred direc- 
tion). This assumption is justified by the fact that the 
results presented in this paper do not change significantly 
if we use p instead of n in the viscosity tensor. 



Throughout our calculations, we will not assume any 
restriction on the viscosity constants except the usual 
req uirements due to thermodynamic stability (see e.g. 
^0| ). Later on we will impose the incompressibility of 
the fluid by assuming a constant mass density p of the 
fluid. We emphasize that this procedure does not require 
any further assumption about the material parameters. 



The set of basic equations is completed by the Gibbs- 
Duhem (the local formulation of the second law of ther- 
modynamics) and the Gibbs relation (which connects the 
pressure P with the other thermodynamic quantities), 
which we will use in the following form: 



de = 



P 



deo + vdg + ipijdVjHi 



h'Arii 



S'(")ds(") + ")dV,s(") 
e + PP + Ta + V ■ g 



(11) 
(12) 



(16) 



—ni + VjVjHi 
1 



(A - l)4nfe + (A + l)4n, V^Vk 



7i 
= 
d 



2 

—^tkhk 



(17) 
(18) 



dt 



+ V^kl^lVk 
-y^P 



Q,(n,s)2(n,s) 



(19) 



(20) 



For the reversible parts of the equations some coupling 
constants have been introduced: The flow-alignment ten- 
sor 



1 



Ayfc = [(A - l)Sljnk + (A + l)4n,] 



(21) 



The newly defined quantities in Eq. (ll|) are connected 
to the thermodynamic forces (Tab. p) by the following 
relations: 



with the flow-alignment parameter A (and using 5,^ 



—riinj) and the coupling between flow and order para- 



meter 



- f3^r^S^^+f3l['^n,n, (22) 
/3j^ - {S,, ~ VrV3 ) + ^i'('v^V3 • (23) 



hi 



h' 



'5e 



ou 



Furthermore there is a reversible coupling between the 

(13) 

layer displacement and the velocity field in equation (p^). 
But its coupling constant has to be unity due to the 

(14) Gallilei invariance of the equations. As mentioned above, 
the use of u is limited to small layer deformations. 

(15) The transverse Kronecker symbols S:lj in Eqs. (17 2^) 
guarantee the normalization of h. This implies that only 
two of the Eqs. (^^ are independent. For the following 
calculations it turned out to be useful to guarantee the 
normalization of the director by introducing two angular 
variables 6 and 6 to describe the director. 



Following the standard procedure within the frame- 
work of irreversible thermodynamics we find the follow- 
ing set of macroscopic hydrodynamic equations [pO| 03, 



sm f cos ( 
sin 9 sin (j 
cos 6 



(24) 
(25) 
(26) 
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Consequently, the Eqs. ( p7| ) have to be replaced using an- 
gular variables. Denoting the right hand side of Eqs. (IT"" 
with Yi, this can be done the following way. 

d_ 

dt 

— Y, sin 



Yt cos 9 cos ( 



Yy cos ( 



d 



Y. 



(27) 
(28) 



In the same way, we guarantee the normalization of p by 
using 



Px 
Py 

Pz 



= 



-VyU 



I -Pi 



(29) 
(30) 

(31) 



The different ways of normalizing fi and p arise from the 
fact, that p is parallel to iz in zeroth order, whereas fi 
encloses a finite angle with for any given shear rate. 

The set of macroscopic hydrodynamic equations we 
now deal with (^ ^ - 1^, |^, H) follows directly from 
the initial input in the energy density and the dissipation 
function without any further assumptions. 

To solve these equations we need suitable boundary 
conditions. In the following we will assume that the 
boundaries have no orienting effect on the director (the 
homeotropic alignment of the director is only due to the 
layering and the coupling between the layer normal p and 
the director n). Any variation of the layer displacement 
must vanish at the boundaries. 



w(±id) = 



(32) 



For the velocity field the situation is a little more com- 
plex: We assume no-slip boundary conditions, i.e. the 
velocity of the fluid and the velocity of the plate are the 
same at the surface of the plate. It is convenient to split 
the velocity field in two parts: the shear field vq which 
satisfies the governing equations and the no-slip bound- 
ary condition and the correction vi to this shear field. 
The boundary condition for vi now reads: 



0. 



(33) 



Making use of the following considerations this condition 
can be simplified. Due to Eq. (Ilq) the z-component of 
Vi is suppressed by a factor of Ap( which is typically ex- 
tremely small jl^, ^ ) . Making use of the results of |2|] 
we can assume that vi depends only on y and z and thus 
conclude [with Eq. (p^j that also the j/-component of vi 
is also suppressed by Ap. For this reason one can assume 
that vi^y and vi^z are negligible and the only relevant 
boundary condition for the velocity field is 



Vl,x = 0. 



(34) 



The validity of this assumption is nicely illustrated by 
our results. Figure ^ shows that Vy and Vz are indeed 
suppressed by Ap. 



TABLE III: If the symmetry under inversion of z is given for 
one component of Xi , the symmetry of all other components 
follows directly from the linearized set of equations. Here we 
give the z-symmetry of all components assuming that u is an 
even function of z. 



Quantity 


z-Symmetry 


Quantity 


2-Symmetry 


u 


even 




even 


e 


odd 


Vy 


odd 




even 


Vz 


even 


p 


odd 




odd 



C. Technique of solution 

The aim now is twofold: Finding a spatially homo- 
geneous solution of the governing equations (for a given 
shear rate) and investigating the stability of this solution. 
In this section we will describe the general procedure and 
give the results in Sec. [II. 

We write the solution as the vector X = 
{9,(l),u,Vrc,Vy,Vz, P, s'^^''^^) consisting of the angular vari- 
ables of the director, the layer displacement, the velocity 
field, the pressure and the modulus of the (nematic or 
smectic) order parameter. For a spatially homogeneous 
situation the equations simplify significantly and the de- 



sired solution Xq can directly be found (see Sec. Ill A). 
To determine the region of stability of Xq we perform a 
linear stability analysis. I.e. we add a small perturbation 
Xi to the homogeneous solutions Xq: X = Xq+Xi (with 
Xi -C Xq) and linearize the governing equations in the 
small perturbations. In short, the solution of the equa- 
tion L Xi = -^^Xi is analyzed. Here L denotes the opera- 
tor for the linearized set of the governing equations. The 
ansatz for the unknown quantities must fulfill the bound- 
ary conditions [see the discussion following Eq. (32)] and 
follow the symmetry scheme given by Tab. Ill, Assuming 



an exponential time dependence and harmonic spatial de- 
pendence of Xi 



Xi^i exp[(ia; - 



T Lsm(gy) 



cos(q^z) 
s\n{qzz) 



(35) 



fulfills all requirements (with an oscillation rate w, a 
growth rate 1/t and a wave vector q = qiy + qz^z)- In 
this ansatz we made use of the results by Oswald and 
Ben- Abraham [ psf , who have shown that in standard di- 
lated smectic A under shear the first instability will set in 
with a wave vector along the neutral direction of the fiow 
{q-Sx — 0). After inserting the above ansatz in the lin- 
earized set of (partial differential) equations, a set of cou- 
pled linear equations is obtained to determine 1/r and lu. 
From the standard smectic A hydrodynamics it is known, 
that shear does not destabilize the layers. Since our ex- 
tended formulation of the smectic A hydrodynamics is 
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equivalent to the standard smectic A hydrodynamics for 
vanishing external fields (e.g. shear rate), we assume 
that the layers are stable for low enough shear rates, i.e. 
l/r < for small shear rates. So 1/r = marks the set of 
external parameters (shear rate) and material parameters 
above which Xi grows. Typically we hold the material 
parameters fixed and the only external parameter is the 
shear rate. The solvability condition of the corresponding 
set of linear equations gives a relation between the shear 
rate [and tilt angle which is directly connected to the 
shear rate, see Eq. ( p8| ) below], uj and the wave vector q. 
For every given q a specific shear rate (and tilt angle 6*0) 
can be determined which separates the stable region (be- 
low) from the unstable region (above). This defines the 
so called curve of marginal stability (or neutral curve) 
Oo{q). If, for any given set of external parameters, the 
tilt angle Oq lies above the curve of marginal stability for 
at least one value of q, the spatial homogeneous state is 
unstable and undulations grow. The smallest shear rate 
(tilt angle) for which undulations can grow is called the 
critical shear rate (tilt angle). Technically speaking, we 
solve L X = ibjX — in many cases we can set uj = Q, see 
below. We point out that this linear analysis is only valid 
at the point where the first instability sets in. Without 
further investigations no prediction of the spatial struc- 
ture of the developing instability can be made. Also the 
nature of the bifurcation (backward or forward) must be 
determined by further investigations. 

For practical reasons we used dimensionless units in 
our numerical calculations. The invariance of the gov- 
erning equations under rescaling time, length and mass 
allows us to choose three parameters in these equations 
to be equal to unity. We will set 

Bi = 1, 71 = 1, and - = = 1 (36) 

and measure all other quantities in the units defined by 
this choice. Nevertheless we will keep these quantities 
explicitly in our analytical work. 

To extract concrete predictions for experimental pa- 
rameters from our calculations is a non-trivial task, be- 
cause neither the energetic constant Bi nor the rotational 
viscosity 71 are used for the hydrodynamic description of 
the smectic A phase (but play an important role in our 
model). Therefore, we here rely on measurements in the 
vicinity of the nematic-smectic A phase transition. Mea- 
surements on low molecular weight liquid crystals made 
by Litster js^ in the vicinity of the nematic-smectic A 
transition indicate that Bi is approximately one order of 
magnitude less than B^. As for 71 we could not find any 
measurements which would allow an estimate of its value 
in the smectic A phase. In the nematic phase 71 increases 
drastically towards the nematic-smectic A transition (see 
e.g. [^). Numerical simulations on a molecular scale are 
also a promising approach to determine these constants 

Due to the complexity of the full set of governing equa- 
tions, we will start our analysis with a minimal set of 



variables (0, (j) and u) a nd supp ress the coupling to the 
other variables (see Sec. [II B 1 ). Step by step the other 
variables will be taken into account. The general picture 
of the instability will turn out to be already present in 
the minimal model, but many interesting details will be 
added throughout the next sections. In comparison to 
our earlier work pS] we now use the way of normaliz- 
ing fi and p derived above. This will lead to some small 
differences in the results but leaves the general picture 
unchanged. First we assume a stationary instability (i.e. 
we let a; = 0); later on we discuss the possibility for an 
oscillatory instability and have a look at some special 
features of the system (Sees. IIIC and HID). 



III. RESULTS AND DISCUSSION 

A. Spatially homogeneous state 

Looking for a spatially homogeneous solution, the gov- 
erning equations simplify significantly. A linear shear 
profile 



(37) 



is a solution to ( [l9|) and u stays unchanged in this regime. 
The only variables which have a zeroth order correction 
for all shear rates are the tilt angle and the modulus of 
the nematic order parameters 5^"^: 

A + 1 



-Asin^(0o) 7 



— — sin(0o) cos((?o) 
71 



^sin(0o)(l-cos(0o)) (38) 
71 

a(")Lo4"^ = (/3j"^-/?l"^)sin(0o)cos(0o)7 (39) 



Equation (^9|) shows that nematic degrees of freedom 
couple to simple shear, but not the smectic degrees of 
freedom; the modulus of the nematic order parameter has 
a non-vanishing spatially homogeneous correction [see 
Eq. (^9|)], whereas the smectic order parameter stays un- 
changed. The reason for this difference lies in the fact 
that /3|"'* and include fi and p, respectively, which 



coupled differently to the flow field [see Eq. ( ^4 g3|)]. 
Eq. (IH) gives a well defined relation between the shear 
rate 7 and the director tilt angle which we will use 
to eliminate 7 from our further calculations. To lowest 
order 9q depends linearly on 7: 



.71 A + 1 
^~b'i 2 



o{el) 



(40) 



We are not aware of any experimental data, which would 
allow a direct comparison with these results. We stress, 
however, that molecular dynamics simulations by Sod- 
demann et al. pQ] are in very good agreement with 
Eqs. (H) and 
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In contrast to the director tilt the lowest order correc- 
tion to the nematic order parameter is quadratic in the 
shear rate (tilt angle). 



(n) 



2 z?!"^-/?!"^ 



O(0 



(41) 



A + 1 71 

In the following we consider perturbations around the 
spatially homogeneous state given above. 



B. Stationary instability 

1. Minimal set of variables 

Let us first consider the effect of our modifications re- 
garding the normalization of n and p in comparison to 
our earlier results . For this purpose we consider only 
a minimal set of variables: the director (characterized by 
the two angles 9 and 0) and the layer displacement u. We 
neglect all couplings of these variables to other quanti- 
ties describing the system, namely the velocity field and 
the moduli of the nematic and smectic order parameters. 
Within these approximations the equations to solve are: 



= 



A 



e |27Asin(6'o)cos(6'o) 



+ ^ \sm^{eo) - cos2(0o) + cos(0o)] 
71 



[sin2(^o)-cos2(^o)] I 
71 J 



Au — siii{9o)qz 
71 



1 B 

Z 7i 



Ae XpBQsm{Oo)qz 

A^ XpBiqsm{6o)cos{eo) 



(42) 



(43) 



An, Xr. 



Bog2(l-cos(0o)) 



Biq cos 



Kq^ + Boqi 



(44) 



Here we inserted an ansatz of the type (35) and denoted 
the linear amplitudes of 9, 4>, and uhy Aq, A^, and 
respectively. One can solve these equations either by ex- 
panding them in a power series of Oq (expecting to get 
a closed result for the critical values) or numerically. It 
turns out, that one has to take into account terms (at 
least) up to order 6*0^ in Eqs. (|4^ - ^) to get physically 
meaningful (but rather long and complicated) analyti- 
cal results. For this reason the closed expressions have 
no advantage over the purely numerical solutions and we 
do not give the analytical approximations explicitly. A 
comparison with the results of Ref. |Q will be given in 
App. We will present and discuss our fin dings us- 
ing the minimal set of variables in Sec. Ill B 2 in direct 



2. Coupling to the velocity field 

In the previous section we have shown that already a 
minimal set of variables supports our picture of the phys- 
ical mechanism. But neglecting the coupling between 
velocity field and nematic director and vice versa is a 
rather crude approximations since it is well known, that 
this cou pling p lays an important roll in nematohydrody- 
namics [|l9[ g^l- t;he natural next step is to include 
this coupling and toperform a linear stability analysis of 



Eqs. 



comparison to the results of the full set of equations. 



. In this case the standard proce- 
dures leads to a system of seven coupled linear differential 
equations. Following the discussion after Eq. (|3|) these 
equations can be solved by an ansatz of the type given in 
Eq. (|35|). This reduces the system of equations to seven 
coupled linear equations which are easily solved using 
standard numerical tools (such as singular value decom- 
position and inverse iteration to find the eigenvectors). 
Due to the complexity of the equations we used Maple 
to determine the final set of linear equations. The key 
ingredients of this Maple script are given in App. |^. 

Figure ^ gives a comparison of typical neutral curves 
for the minimal model and calculations including the ve- 
locity field. The overall shape of the neutral curve is not 
changed due to the coupling to the velocity field but a 
shift of the critical values (especially in the critical tilt 
angle) is already visible. The inset shows the relative 
amplitudes of the linear solutions at onset on a logarith- 
mic scale. For 9, </> and u the left bars correspond to 
the minimal model and the right bars to the extended 
version. Note that amplitudes with a different sign are 
shown with a different line style in the histograms (see 
figure caption for details). 

Let us have a closer look at the differences between 
the minimal and the extended set of equations and fol- 
low these differences along some paths in the parameter 
space. As mentioned in Sec. [I C, we can omit some of the 
physical parameters by using dimensionless parameters. 
In Figs. ^ ~ ^ we show the dependence of the critical val- 
ues of the tilt angle and wave vector on the dimensionless 
parameters [as defined in Eq. (|36|)]. For all these figures 
we used the same basic set of parameters: Bq = 10, 
K = 10~^, A = 1.1, 1^1 — 1^2 = >^3 ^ 1^4 — = 0.1 and 
Xp = 10~^. These values are estimates for a typical ther- 
motropic LMW liquid crystal, where we made use of the 
results of |3^, |3^ (as far as Bi and 71 are concerned, see 
also the last paragraph in Sec. IIIBl). For flow align- 
ment parameters in the range 1 <; A < 3 the critical 
values vary strongly with A (see Fig. Therefore we 
discuss in addition the situation for A = 2 to indicate the 
range of possible values. 

Considering the critical values as a function of the com- 
pression modulus Bq results in a rather simple situation 
(Fig. 1^): For small values of Bq a significant infiuence of 
the coupling between the director and velocity field is ap- 
parent, which also shows a strong dependence on A. For 
large Bq all these differences vanish and only one single 
curve is obtained. At this point a comparison to dilated 
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FIG. 4: A typical picture for the comparison of the neutral 
curves using the minimal set of variables ( ) and includ- 
ing the velocities field ( ) . The overall behavior does not 

change but the critical values are altered due to the coupling 
with the velocity field. For this plot we used (in the dimen- 
sionless units discussed in Sec. |lCJ) Bo = W, K — 10^^, 
X = 1.1, ui = Vi = 1^3 = ^4 = 1^5 = 0.1 and Xp = 10~®. 
The inset shows the linear amplitudes Ai (where i stands for 
9, 4>, etc) at onset. Since the logarithm of the amplitudes is 
shown, amplitudes with different sign are shown with a differ- 
ent line style. Using the minimal set (left bars) all amplitudes 

have the same sign ( ). Including the velocity field (right 

bars) some amplitudes are positive ( ) , others negative 

( ) . Note that we use in this and all following plots the 

dimensionless units defined by Eq. (Bq) . 




FIG. 5: A significant difference between the various ap- 
proaches is only visible for Bo < 100. At higher values of 
Bq the number of free variables plays no noticeable role and 
the critical values follow a master curve. The solid lines 

( ) show results including the velocity field, the dashed 

lines ( ) correspond to the minimal set of variables. At 

low Bo in the upper curves we used A = 2. 
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smectic A is instructive. It is well known [|T, 22 that in 
dilated smectic the critical wave vector and the critical 
dilatation show a power law behavior as a function of 
with exponents 1/4 and —1/2, respectively. In the limit 
of large -Bq we found the same exponents already in our 
earlier analysis p5|. If we fit power laws to our results 



FIG. 6: Plotting the critical values as a function of the bend- 
ing modulus K shows a convergence of the curves, which is 
nevertheless not as pronounced as in the case of Fig. |. The 
influence of A on the critical tilt angle is significant (A = 2 
in the upper curves and A = 1.1 in the lower ones). Again 

the solid lines ( ) show results including the velocity field 

and the dashed lines ( ) correspond to the minimal set of 

variables. 
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FIG. 7: In all our calculations vi^x is the dominating compo- 
nent of v\. This graph demonstrates that the other compo- 
nents are suppressed by Ap (making them almost negligible). 




for i?o > 10^ we find the exponents equal to « 0.235 
and K, —0.37 for qc and 0c, respectively (note that the 
dilatation in our model is « ^6*^). So both approaches 
(the minimal set of variables and the calculations includ- 
ing the velocity field) show, despite all similarities to the 
standard model of smectic A and to our earlier analysis, 
differences in the details of the instability. 

A similar, but less pronounced, situation is apparent, 
when plotting the critical values as a function of the 
bending modulus (see Fig. ^. The curves tend to con- 
verge for large K , but there remains a difference between 
the minimal set of variables and the calculations includ- 
ing the the velocity field. Fitting the if-dependence with 
power laws (here for K > 10^"*) only the critical wave 
number exhibits an exponent close to the values expected 
from dilated smectic A (« —0.26 vs. —\)- This illus- 
trates the fact that shearing a lamellar system is similar 
to dilating it but not equivalent. 

In contrast to the cases discussed above, the perme- 
ation constant Xp has no strong influence on the critical 
values. For dimensionless values A„ < 10~^ the critical 
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values do not change at all with Ap. For large values 
variations within a factor of two are possible. The per- 
meation constant is known to be very small. In our di- 
mensionless units we expect it be of the order of < 10^^ 
for LMW thermotropic liquid crystals and neglect its in- 
fluence on the critical values for this reason. In Sec. [I B 



we have emphasized that the y- and z-components of the 
velocity field are suppressed via Ap. These qualitative ar- 
guments are clearly confirmed by our numerical results: 
In all our calculations Vi^x is the dominating component 
of Vi and the ratio is of the order of Ap over 

the whole range of physical relevant values of Ap (see 
Fig. This fact nicely supports our argument that we 
can neglect the boundary condition for vi.y^ because vi^y 
vanishes anyway. 

Out of the five viscosities only two (1^2 and v^) show 
a significant infiuence on the critical values. In Fig. ^ 
we present the dependence of 9c and qc on an assumed 
isotropic viscosity (upper row) and on these two viscos- 
ity coefficients (middle and lower row). Since the flow 
alignment parameter A has remarkable influence on these 
curves we have chosen four different values of A in this 
figure, namely A = 0.7, A = 1.1, A = 2 and A = 3.5. 
The curves for A < 1 and A > 3 for an isotropic viscosity 
tensor are very similar to the corresponding curves where 
only 1/2 is varied. In this parameter range the coefficient 
V2 dominates the behavior. Note that the influence of 
on the critical values is already much smaller than that 
of V2. We left out the equivalent graphs for the other vis- 
cosity coefficients, because they have almost no effect on 
the critical values. For further comments on the influence 



of an anisotropic viscosity tensor see also Sect. |IID. 

All the parameters we have discussed up to now caused 
variations in the critical values that did not select speciflc 
values of the considered parameter. In this aspect the 
situation is completely different in the case of the flow 
alignment parameter A. As shown in Fig. ^ there is a clear 
change in behavior for A « 1 and A « 3. The critical tilt 
angle is increased for values of A in this interval and the 
critical wave vector tends to rise only at the boundaries of 
the interval and is reduced in between. Fig. ^ illustrates 
how this structure depends on the viscosities (assuming 
all five viscosities to be equal) and on the elastic constants 
of the layers. In the first row we follow this behavior 
for viscosities varying from Vi = 1 down to Vi — 10^^. 
Clearly, the influence of A is more pronounced the lower 
the viscosities are. Both elastic constants of the layers, 
the compressibility Bq and the bending modulus K (in 
our dimensionless units Bi — 1), have in general a similar 
influence on the shape of the graphs: the smaller the 
elastic constants are, the more pronounced the structure 
becomes. For this reason we just give the plot for So 
(second row in Fig. ^) and omit the plot for K. 

These dependencies on the system parameters give 
some important hints for an interpretation of Fig. ^. The 
currents and quasi-currents for the velocity field and the 
director consist of two parts [see Eqs. ( p7| ) and (|Tg|)]: a 
diagonal one (coupling e.g. the components of v among 



FIG. 8: Only the viscosities V2 and 1/3 can influence the critical 
parameters significantly. The upper row depicts the depen- 
dence on a isotropic variation of the viscosity. In the middle 
and lower row we present the variation with V2 and 1/3 setting 
the other viscosities to Ui = 0.1. Here the thick solid lines 

( ) represent the minimal set of variables. For the full 

set of variables we have chosen four different values of A: the 

solid curves ( ) with A = 0.7, the dashed curves ( ) 

with A = 1.1, the dotted curves ( ) with A = 2 and the 

dot-dashed curves ( ) with A = 3.5. Note the similarities 

between the curves for small ( ) and large A ( ) in the 

upper and middle row: In these regimes U2 is the dominating 
viscosity. 
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each other) and an ofl-diagonal one (coupling the direc- 
tor to the velocity held). The former ones are propor- 
tional to the elastic constants or to the viscosity tensor 
whereas the latter one is a function of the flow alignment 
parameter. So reducing either the elastic constants or 
the viscosities increases the portion of the cross-coupling 
terms in theses equations. I.e. the observed tendencies 
are exactly what one would expect. The next step in the 
interpretation of the shape of the curves is to have a closer 
look at the structure of the cross-coupling term. The flow 
alignment tensor X^jk = 5 [(A — 1)51^11^. -I- (A -I- l)i5,^nj] 
obviously changes its behavior for A = 1: The first part 
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FIG. 9: Plotting the critical values as functions of the flow 
alignment parameter A reveals an interesting structure for 
1 ^ ^ 3. In the upper row we plot this dependence for 
a set of (isotropic) viscosities ranging from Ui = 1 (thick solid 

line, ) down to Vi = 10"'^ (thick dashed line, ). 

The dotted line ( ) reveals that this dependence on A is 

absent in the minimal model. The lower row illustrates the 
behavior for varying layer compressibility Bo with Bo ~ 3 for 

the thick solid curve ( ) and Bq = 100 for the thick dashed 

curve ( ) . In all plots the thin solid lines ( ) give the 

behavior for some intermediate values. For an interpretation 
of this behavior see the text. 





assumed that varies only slightly around Sq"'^-* and 

thus only the lowest order terms in s*^"'*' contribute to 
the free energy. For the spatially homogeneous state we 
had [see Eqs. (^ |4l| )] a correction to the nematic S*'"' 
proportional to the square of the shear rate (flg 7 for 
low 7): 



A + 1 



71 



(") 



0(^0) 



As a consequence s^f^ must be small compared to 5q"'*^ 
(which is by construction limited to the range < 
g{n,s) Thus a reasonable restriction is 



c(")i 



< 0.5. 



(45) 



As shown in Fig. [l^, evaluating this relation at the on- 
set of the instability reduces significantly the physically 
reasonable range for some parameters. This restriction 
applies only for the nematic material parameters and, 
in general, nothing can be said about the corresponding 
smectic parameters. We will, however, take the smectic 
parameters in the same range as the nematic ones. If not 

= 0.1 
10-4, 



indicated otherwise we used Lq"'"-* 



/3| 



(ti,s) 
I 



0.005, 



Mo = 

Q,(",^) ^ 



, 4"''^ = 0.01, 
/3^"'^^ = 0.01, 



0.001 for the plots of 



this section (along with parameter set specified in the 
previous section). 



The ansatz for s\ 



(n,s) 



following Eq. 
1 



5^ reads 



A*"'^^ exp[(iw + -)t] sin(q^z) cos{qy). (46) 



changes its sign. Note that we are in a region of the 
parameter space, where Xijk is a dominating term (since 
either the viscosities or the elastic constants are small). 
Additionally, 6:^nk contains up to the third power of one 
director component. For this reason we expect that — 
in the linearized set of equations — some coupling terms 
change their sign for A = 1 others for A = 3. E.g. the (p- 
component of the director is coupled to the x- and z com- 
ponent of the velocity field by the terms (A — l)/2 dyVx 
and (A— 1)/2 cot{do)dyVz- Similarly the reversible part of 
the coupling of Vy to vanishes for A = 3. The monitored 
structure in the plots cannot be attributed to one single 
cross-coupling term, but the given examples demonstrate 
that something should happen in this parameter range. 



3. Including the order parameters 

In the preceeding paragraphs we investigated undula- 
tions assuming a constant modulus of the order parame- 
ter 5*^"'*^ = 5'o"''*'' -|- s^P'^K In general one would expect 
that the undulations in the other observable quantities 
should couple to some extent to the order parameter. In 
the formulation of the free energy (see Sec. II B) we have 



The modulations of 5*^"''*-' in the linear analysis are max- 
imum at the boundaries and in phase with the layer dis- 
placement u. The sign of the amplitude ^i"'^' depends 
on the coupling to the velocity field (only the anisotropic 

part — Z?^-* is relevant) and on the coupling to the di- 
rector undulations (via Mijk, only for the nematic ampli- 
tudes aI"^). If one assumes that shear reduces (and does 
not increase) the modulus of the order parameter, the 



nematic /3||"'' — I3'"J^' is positive [Eqs. ( |39D and ([4l|)]; once 

(s) (s) 

again nothing can be said about the smectic — . 

In general the critical values are not at all or only very 
slightly infiuenced by the coupling to the modulus of the 
order parameter (see Fig. pd] ). Figure |ll| summarizes the 
parameters with the largest infiuence on ^i"'*^ . In almost 
all investigated cases the modulation of the nematic or- 
der is much larger than in the smectic order. Whether 
the order is reduced or increased in regions where the 
layers are compressed depends in the phenomenological 
constants /Jy"'*"* — /j'j"'*'' and Mq which have not been 
measured up to now. 

The above results reveal some interesting features. As 
shown in Tab. Ill, the modulations of the order parame- 



ter change sign under inversion of the z-axis. Considering 
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FIG. 10: Evaluating Eqs. @ and at onset gives an im- 
portant restriction on the range of possible parameter values 
(here the cases of a'"' and /Jy"' — /3*j"'). Note that the critical 
9o is a function of the material parameters. 
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FIG. 11: Out of the material parameters connected with the 
order parameter, only p'^^"'"^ — p'^^'"'^ has a measurable effect 
on the critical values. Some more parameters can influence 
the amplitudes of the order parameter undulation, namely 
Z/'i"'^' and Mo (the latter one is only present in the case of 
the nematic order parameter) . All amplitudes have been nor- 
malized such that = 1. Note that the smectic Ai"^ has 
been multiplied by 10^ in the right column. For a better 
comparison we used a log-log scale in the lower left plot and 
changed the sign of jA^"' in this plot. 
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the boundary condition (i.e. taking our ansatz) this leads 
to the fact that the effect on the modulus of the order 
parameters is maximum at the boundaries. So the linear 
analysis predicts that the regions where the order para- 
meter is influenced most by the undulations are close to 
the boundaries. Since the probability for the formation 



FIG. 12: In most parts of the scanned parameter space no 
possibility for an oscillatory instability was found. If the di- 
rector field is only very weakly coupled to the layering (in 
this plot we used Bo = 50 and Vi = 5) a. neutral curve for 
an oscillatory instability ( ) appears above the station- 
ary neutral curve ( ). Note that the critical wave vectors 

are close to each other for both, oscillatory and stationary 
instability. The inset shows the frequency along the neutral 
curve. 




of defects is higher in places where the order parameter 
is lower, we have identified areas where the creation of 
defects is facilitated. But our analysis does not allow 
to predict the structure of the defects. Nevertheless this 
effect gives a possible way how to reorient the parallel 
layers. Interestingly, experiments in block copolymers 
by Laurer et al. ||^ show a defect structure close to the 
boundaries which is consistent with this picture. 



C. Oscillatory instability 

All our arguments in the previous sections were based 
on the assumption that the undulations set in as a sta- 
tionary instability. I.e. that the oscillation rate uj in our 
ansatz Eq. (^5|) vanishes at onset. In this section we will 
discuss the situation for non-zero uj and find that our 
previous assumption was justified. In our linear analy- 
sis enters now (for the first time in this paper) the mass 
density of the system, which we will choose to be equal 
to unity p — \. 

The search for a possible oscillatory instability is 
slightly different from the procedure used in the station- 
ary case. The solvability condition of the linearized set 
of equations determines both the neutral curve and the 
frequency along this curve. When searching for such a 
solution we scanned approximately the same parameter 
space as used for Figs. Since the frequency tends to 

zero when the oscillatory neutral curve gets close to the 
stationary one, we concentrated on the frequency range 
< a; < 2 and check in some cases for higher frequencies. 
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It turned out that only in cases when the director field 
is very weakly coupled to the layering a neutral curve for 
an oscillatory instability is possible. This weak coupling 
manifests itself in small Bi and 71, which is in our set 
of dimensionless variables equivalent to large Bq and Vi. 
Oscillatory neutral curves were only found for Bq > 100 
or Vi > 1. In all investigated cases a oscillatory neutral 
curve is either absent or lies above the neutral curve for 
a stationary instability. When a oscillatory neutral curve 
is possible, it ends in the points where it meets the sta- 
tionary neutral curve (see Fig. |l2|). The corresponding 
frequency approaches zero in the end points of the os- 
cillatory neutral curve. If we ignore for the moment the 
stationary neutral curve and consider only the oscillatory 
instability, the corresponding critical values are found to 
be rather close to the stationary one and to approach 
them the weaker the coupling between the director and 
the layers becomes. To summarize, an oscillatory insta- 
bility was not found to be possible in all investigated 
cases and seems to be extremely unlikely to occur. 



D. Anisotropic viscosity 

In Fig. H we have illustrated that a small viscosity co- 
efficient 1/2 facilitates the onset of undulations. In this 
section we will have a closer look at the effect of an 
anisotropic viscosity tensor and ask whether undulations 
can be caused only due to viscosity effects without any 
coupling to the director field (i.e. we consider standard 
smectic A hydrodynamics in this section). 

Let us start our considerations by looking at the spa- 
tially homogeneous state. In a sample with parallel align- 
ment the apparent viscosity is 1/3, which can easily be 
seen from the force on the upper boundary: 



I = 62 • CT = 71/361; 



(47) 



Similarly the viscosity of a perpendicular alignment is 
given by h'2- 



Fi = e 



z ■ g_ = 71/263; 



(48) 



For 1/2 < 1/3 a simple shear flow in a perpendicular align- 
ment causes less dissipation than in a parallel alignment. 
The next step is to study the stability of these alignments 
in the linear regime. Following the standard procedure 
(as described above) we find a solvability condition of the 
linearized equations which does not depend on the shear 
rate 7. 



= 



[Boqi 



V3 {q^ -qlf + '2 (t^2 + J^s) q^ql] } X 
Kq^) (1/29' + i^sql) (49) 



Consequently, a parallel alignment of smectic layers is lin- 
early stable against undulations even if the perpendicular 
alignment might be more preferable due to some thermo- 
dynamic considerations. As we have shown in Fig. |^, this 



rigorous result of standard smectic A hydrodynamics is 
weakened in our extended formulation of smectic A hy- 
drodynamics. When the director can show independent 
dynamics, an appropriate anisotropy of the viscosity ten- 
sor can indeed reduce the threshold values of an undula- 
tion instability. 



IV. COMPARISON TO EXPERIMENTS AND 
SIMULATIONS 

In the previous sections we have shown that the in- 
clusion of the director of the underlying nematic oder in 
the description of a smectic A like system leads to some 
important new features. In general, the behavior of the 
director under external fields differs from the behavior of 
the layer normal. In this paper we only discussed the ef- 
fect of a velocity gradient, but the effects presented here 
seem to be of a more general nature and can also be ap- 
plied to other fields. The key results of our theoretical 
treatment are a tilt of the director, which is proportional 
to the shear rate, and an undulation instability which 
sets in above a threshold value of the tilt angle (or equiv- 
alently the shear rate). 

Both predictions are in agreement with experimental 
observations. For side-chain liquid crystalline polymers 
Noirez [|l^ observed a shear dependence of the layer 
thickness. In the parallel orientation the layer thickness is 
reduced by several percent with increasing shear. To our 
knowledge, two groups have investigated the evolution 
of a parallel alignment to the vesicle state for lyotropic 
systems (see Miiller et al. Q and Zipfel et al. Q). In 
both papers the authors argue that cylindrical structures 
(with an axis along the flow direction) are observable as 
intermediates. These observed cylindrical intermediates 
are very close to the undulations proposed by our theo- 
retical treatment. 

For an approximate quantitative comparison of our 
theoretical results with the experiments on lyotropic liq- 
uid crystals we make a number of assumptions about the 
material parameters. As we have shown in Sec. [II B the 



different approaches cause only small variations in the 
critical wave number. For this estimate it suffices to use 
the c ritic al wave number obtained in our earlier work [see 
Eq. ([A^)]. For lyotropics it is known that the 

elastic constants can be expressed as 



and 



B = ^n^ — 
64 K 



{1-5Y 



(50) 



(51) 



where k = a^fesT is the bending modulus of a single 
bilayer, I is the repeat distance, 5 is the membrane thick- 
ness, kB is the Boltzmann constant, T is temperature and 
a„ is a dimensionless number of order of unity. With this 
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relations we can estimate the cr itic al wave vector for a 
sample of thickness d using Eq. ([A^): 



I 



3^2 



8a^d {l-6y 



(52) 



The parameters of the experiment by Zipfel et al. ||4l|] 
are: d = 1 mm, S = 2.65 nm, I = 6.3 nm and = 
1.8 On this basis we estimate the critical wave 

length to be of the order of 



A. 



6.4 /im 



(53) 



Zipfel et al. |4l|] observed a vesicle radius of 3 /im, which 
is clearly compatible with our calculation. We note that 
this estimate assumes that the experiments are done in 
the hydro dynami c regime. 

In Sect. IIIB 3 we have pointed out that the effect on 
the order parameter is maximum close to the boundaries 
of the layer. In a reoriented sample Laurer et al. ||] 
have identified defects near the boundary of the sample 
which are in accordance with the predicted influence on 
the order parameter. 

Molecular dynamic simulations recently made by Sod- 
demann et al. psf offer a very precise insight in the 
behavior of the layered systems under shear. Direct com- 
parison of these simulations to the analytic theory pre- 
sented above show a very good agreement between both 
approaches pOf . 

The mechanism we have proposed here is somewhat 
similar to a shear induced smectic-C like situation. Con- 
sequently, undulations should also be observed near the 
smectic-A-smectic-C transitions. Indeed, Johnson and 
Saupe ||46|| and later Kumar ||47|| report such undulations 
just below the transition temperature. In the same spirit 
Ribotta and Durand Q report a compression induced 
smectic-C like situation. 

To conclude, we have shown in this work that the in- 
clusion of nematic degrees of freedom in the description 
of smectic A like systems opens the way for a shear in- 
duced destabilization of the layers under shear. Our re- 
sult are compatible with experimental observations and 
are in good agreement with molecular dynamics simula- 
tions. 
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APPENDIX A: MINIMAL ANALYTIC MODEL 

In our earlier work [p5| we considered two independent 
variables (the layer displacement and the y-component of 



the director) and found the critical values to be 

' and 




(Al) 
(A2) 



To compare our present analysis to these results we ex- 
pand Eqs. (^, ^4|) in power series in Oq (up to 9q) and 
take only the terms connected with and u. 



Ad, — ^^'o 

71 



AiiXp 



Au—q 
71 



(A3) 



Biq^ + Kq^ + Boql - ^O^q^iBo + 2B,] 



A^XpBiOoq 



(A4) 



The solvability condition of Eqs. (A3, A4) defines the 
neutral curve Oo(q) and its minimum directly gives the 
critical values 9c and qc (within the approximations of 
this section). 

(A5) 
(A6) 




The differences between Eqs. (Al, A2) and (A5 - A7) are 
mainly due to the correct normalization of p [see Eqs. (^ 
- 13^)] used in the present paper. 

To summarize, we conclude that our former results are 
a special case of the present analysis when the correct 
normalization of p is implemented. Especially the diver- 
gence of the critical values at Bq — 2Bi turns out to be 
an artifact of the normalization of p used in Ref. p5| . 



APPENDIX B: GENERATING THE SET OF 
LINEAR EQUATIONS 

Since the theoretical methods used in this paper (irre- 
versible thermodynamics and linear stability analysis) of- 
fer well defined algorithms for the generation and analysis 
of macroscopic hydrodynamic equations, we performed 
parts of the calculation using Maple. In this appendix 
we describe the key ingredients of a suitable Maple pro- 
gram. A good starting point for such an approach are 
the balance equations for the unknown quantities ( |l6| , 
|l| - 1^, |2^, and H) along with the energy density ^) 
in the appropriate approximation. These equations are 
entered directly in Maple, with the unknown quantities 
being functions of time and the spatial coordinates. The 
thermodynamic forces used in these equations are deter- 
mined by Eqs. (13 -|l^). For an implementation of these 



equations one must take into account that Maple can only 
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compute the derivative with respect to constants and not 
with respect to functions, i.e. the relevant functions in 
the energy density must be substituted temporarily by 
constants. 

For the linearized set of equations we substitute the 
unknown quantities by expressions of the type 

6{t, X, y, z) = 6q + aAg sm{qzz) cos{qy) exp{iu}t) (Bl) 

in the governing equations. Here a is a small parameter 
and Ae is the relative amplitude of the linear correction 
to ^0- Expanding the substituted set of equations in a 



power series in a gives in zcroth order the spatially ho- 
mogeneous equations and in first order the linear set of 
equations, which arc no longer differential equations but 
algebraic ones. One obtains a matrix representation of 
these equations by expanding them in power series of the 
relative amplitudes and taking only the first order terms. 
After dividing by the terms which depend on the spa- 
tial and temporal coordinates the Fortran code of this 
matrix representation is generated using the CODEGEN, 
FORTRAN function of Maple and subsequently solved us- 
ing standard numerical procedures. 
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